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Abstract 

ti ■ A voting model (or a generalization of the Glauber model at zero temper- 

ature) on a multidimensional lattice is denned as a system composed of a 
C^ . lattice each site of which is either empty or occupied by a single particle. 

The reactions of the system are such that two adjacent sites, one empty 

the other occupied, may evolve to a state where both of these sites are 

' O ■ either empty or occupied. The continuum version of this model in a D- 

dimensional region with boundary is studied, and two general behaviors of 
such systems are investigated. The stationary behavior of the system, and 
the dominant way of the relaxation of the system toward its stationary 
state. Based on the first behavior, the static phase transition (discontin- 
£NJ , uous changes in the stationary profiles of the system) is studied. Based 

^ ' on the second behavior, the dynamical phase transition (discontinuous 

C*) ' changes in the relaxation-times of the system) is studied. It is shown that 

C/^ . the static phase transition is induced by the bulk reactions only, while the 

' S SL ' dynamical phase transition is a result of both bulk reactions and boundary 

25^ , conditions. 
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1 Introduction 

The study of the reaction-diffusion systems, has an attractive area. A reaction- 
diffusion system consists of a collection of particles (of one or several species) 
moving and interacting with each other with specific probabilities (or rates in the 
case of continuous time variable). In the so called exclusion processes, any site 
of the lattice the particles move on, is either vacant or occupied by one particle. 
The aim of studying such systems, is of course to calculate the time evolution 
of such systems. But to find the complete time evolution of a reaction-diffusion 
system, is generally a very difficult (if not impossible) task. 

Various methods have been used to study the reaction-diffusion system: an- 
alytical techniques, approximation methods, and simulation. The success of the 
approximation methods, may be different in different dimensions, as for exam- 
ple the mean field techniques, working good for high dimensions, generally do 
not give correct results for low dimensional systems. A large fraction of analyt- 
ical studies, belong to low-dimensional (specially one-dimensional) systems, as 
solving low-dimensional systems should in principle be easier [1—12]. 

Various classes of reaction-diffusion systems are called exactly-solvable, in 
different senses. In [13-15], integrability means that the V-particle conditional 
probabilities' S-matrix is factorized into a product of 2-particlc S-matrices. This 
is related to the fact that for systems solvable in this sense, there are a large 
number of conserved quantities. In [16-24], solvability means closedness of the 
evolution equation of the empty intervals (or their generalization). 

Consider a reaction-diffusion system (on a lattice) with open boundaries. 
By open boundaries, it is meant that in addition to the reactions in the bulk of 
the lattice, particles at the boundaries do react with some external source. A 
question is to find the possible phase transitions of the system. By phase tran- 
sition, it is meant a discontinuity in some behavior of the system with respect 
to its parameters. Such discontinuities, may arise in two general categories: in 
the stationary (large time) profiles of the system, and in the time constants de- 
termining the evolution of the system. In the first case, static phase transitions 
are dealt with; in the second case, dynamical phase transitions. For a review on 
dynamical phase transitions, one can see for example [25]. 

There are systems for them the equation of motion of the one-point function 
(the probability that a certain site be occupied) is closed, that is independent of 
the more-point functions [26-28]. Among these systems is the so called voting 
model (or a generalization of the Glauber model at zero temperature) . In [29] 
a voting system on a one-dimensional lattice was studied, for which at the 
boundaries of the lattice there are injection or extraction of the particles. Based 
on the evolution of the one-point functions, it was shown there that the system 
exhibits two kinds of phase transitions: a static phase transition, corresponding 
to a discontinuous change in the stationary profile of the one-point function; 
and a dynamical one, corresponding to a discontinuous change in the behavior 
of the relaxation time of the system toward its stationary state. In [30-34] , the 
phase structures of extensions of such systems on a one-dimensional lattice were 
investigated. All of these are restricted to the case of a one-dimensional lattice. 



Extending these investigations to higher-dimensional cases would be inter- 
esting. Here, we want to study a multi-dimensional extension of the voting 
model, on a continuum rather than a lattice. 

The scheme of the present article is as follows. In section 2, the multidimen- 
sional voting model on continuum is presented, and the evolution equation for 
the density of the particles is obtained. In section 3, the time-independent of the 
system is studied and it is shown that system exhibits a static phase transition. 
In section 4, the relaxation of the system toward its stationary state is studied 
and it is shown that the system exhibits a dynamical phase transition. Section 
5 is devoted to the concluding remarks. 

2 Multidimensional voting models on continua 

In [11] and [29], a one dimensional voting model (or a generalization of the 
Glauber model at zero temperature) on a lattice was defined as follows. Let 
the system consist of a one-dimensional lattice, each of the sites of which are 
either empty (0) or containing a single particle (A), and let there be a reaction 
between two neighboring sites like 

A0 — > A A, with the rate u + , 

$A -> 00, with the rate u+, 

%A — > AA, with the rate w_, 

A® -> 00, with the rate u— (1) 

In [11], an open lattice was investigated while in [29], a lattice was studied at 
the boundaries of which injection and extraction of particles could occur. It was 
shown that these models are autonomous, meaning that the evolution equation 
of the n-point functions contain only n- or less-point functions. In fact, as it 
was seen in [29], 

— (m) = u + (n 4 _i) + u- (rii+i) - (u+ + u-) («,-). (2) 

Here m is the particle number operator at the site i of the lattice. 

Now consider a multi-dimensional lattice, each site of which is either empty 
or occupied by a single particle, and let there be a reaction like 

A0 — ► AA, with the rate itj, 

%A -> 00, with the rate u h . (3) 

Here, we are considering the reaction between a site i (the right site), which 
is the ending point of the link I, and another site (the left site) which is the 
starting point of the same link. In a one-dimensional lattice, each site is the 
ending point of two links, which had been denoted by + and — . From lp|l. it is 
seen that the evolution equation for the one-point function is 

-jT ( n t) = ^i u l (( n i-l) ~ ( n t-l n *)) _ u l (( n i) - ( n i-l n i))}> ( 4 ) 



where by the site index i — I, it is meant a site which is the starting point of 
the link I, the ending point of which is the site i. It is seen that the two-point 
function in the right-hand side of 10} cancel each other. So, 

-jl (ni) = J2ui ({rii-i) - (ru)). (5) 

i 

Now, assume that the one-point function is a slowly-varying function of its 
argument (i). In this case, one can define a smooth particle density function of 
the continuous position variable r, with 

p(ri):=-{m}, (6) 

where i*j is position of the lattice site i, and V is the specific hypervolume of a 
site. Then (JjjJ can be rewritten as 



d_ 
di 



I 



Si ■ V + \ (Si ■ V) 2 



P, (7) 



where Si is the link-vector, equal to the position of the ending point of the link 
I minus the position of the starting point of the link I, and higher-derivative 
terms have been neglected. Using suitable coordinates for r, one can write the 
second derivative as 

i£»,(VV)* = £(A) 2 , (8) 

where x a: s are the coordinates of r. So, Q) is rewritten as 

|p=(-v-V+V>, (9) 

where 

v:=^u;<Jj. (10) 

i 

Eq. (O is nothing but a diffusion equation combined with a drift velocity v. 

Suppose that © holds for the interior of the region V. Integrating @ on 
V, one arrives at 

— fdVp=-(f dSn-vp+(f dSn-Vp. (11) 

d£ Jv JdV JdV 

The first term in the right-hand side is the rate of change of the total number 
of the particles inside, as a consequence of the drift, while the second term is 
the effect of injecting or extracting particle at the boundary. The boundary 
condition 

n • V/o = a — (3 p, at the boundary (12) 



corresponds to an injection rate of a per unit hyperarea of the boundary, and 
an extraction rate of /3 per unit hyperarea per particle density at the boundary. 
In general, one can take a and (3 position-dependent. 

Comparing © and (|12[1 with eq. (7) of [29], it is seen that one can transform 
eq. (7) of [29] to @ and JHJ) through 



5= .. 

u + v 

v -x = 6 (u — v) 



(3- = 



SuV 

a + a' 

Su 
b 



dvV 



f>+ - >-£, (13) 

V 

where the right-hand sides are the quantities defined in [29] , and the superscripts 
— and + refer to the left- and right-boundaries, respectively. 

From now on, for simplicity we restrict ourselves to the case that the volume 
V is a D-dimensional hyperball with radius R, the boundary of which is a 
hypersphere. 



3 The time-independent state and the static phase 
transition 

Let po be the time- independent solution to © and (|12|l . Using the ansatz 

F q (r) := exp(q • r) (14) 

(with q a constant vector) as a time- independent solution to El one arrives at 

qq-vq = 0, (15) 

which leads to 

q=i(v + v'), (16) 

where v' is an arbitrary constant vector subject to the condition 

(17) 



So, one can write the general time-independent solution to © as 

p Q (r)=Jdn'A(n')F (l (r) 7 

= ydO'A(0')exp|i[(v + v')-r-|v + v'|ii]|, 

=: f dO' A(n') exp[G(v', r)], (18) 

where il' denotes the angular coordinates of v', and A is an arbitrary function. 
It is easy to see that the maximum value of G is zero, and this maximum value 
is reached at a point on the boundary (r = R), where r is parallel with v + v'. 
For large values of R and r, G is a rapidly-varying function and the integral 
is mainly determined from that point of the integration region which maximizes 
G. Generally, there may be two such points: One point is 

(19) 

The other point is 

(v 2 +v) -r = |v 2 +v|r, for r = R, (20) 

which means that q is parallel with r. As the angle between q and v cannot 
exceed 7r/2, the second point exists only if the angle between r and v is less 
than 7r/2. One has 

G(vi,r) = 0, (21) 

and 

G[v 2 (r),r]=G[v 2 (i?r/r),r]+0[(i?-r) 2 ], 

0[(r-Rf. 



_(r-i?)|v + v' 2 | , „ u _ ,. 



2 
(r — R) v • r 



R 

Using (|21|l and 122|) . one arrives at 



0[(r-R) 2 }. (22) 



Ci(fi), r-i?,r-v<0, 

Ci(n) + C 2 (n)exp[ (r -g v - r ], r^i?,r-v>0, 



From this, 

Vp (r = i?) cxn(n-v)6»(n-v), i? -> oo, (24) 

Where is the step function. It is seen that in the thermodynamic limit (i? — + 
oo), the density profile at the boundary is stationary, unless v • r > 0. So, 
changing v one can induce a discontinuous change in the slope of the profile 
density at the boundary. This is the static phase transition, which is seen to be 
independent of the injection and extraction terms, but dependent on the drift 
velocity. 



4 The relaxation of the system toward the sta- 
tionary state, and the dynamic phase transi- 
tion 

Starting from @ and JT2J, one arrives at 
d 



01 



(p-po) = (-v-V+ V z )(p- Po ), 
= : h(p-po), 



and 



n • V(p — po) = —j3 (p — po), at the boundary 



(25) 



(26) 



where po is the time- independent solution to JjJJ and l|12|) . Let ^ be an eigen- 
function of h corresponding to the eigenvalue E. Using the ansatz l|14|l in the 
eigenvalue equation corresponding to h, one arrives at 



which leads to 



qq-vq = £!, 



1 / 
q= _( v + v ') ; 



(27) 
(28) 



where v' is an arbitrary constant vector subject to the condition 

v'-v' = v-v + 4£. (29) 

So, one has 



V>(r) = exp(v • r/2) f dfl' A(Q') exp(v' • r/2), 



(30) 



where A is to be found so that the boundary condition (|26|l is satisfied with ip. 
If the righthand side of l|29|) is positive, then v' is real and for large r, one 
can approximate %j) like 



i[)(r) - exp(v • r/2) A(O) exp(w' r/2), 



(31) 



where i7 is the angular coordinates corresponding to r. The boundary condition 
H16(l . then becomes 



p(Q) 



A(Q) = 0. 



(32) 



This has a nonzero solution for A, provided the parenthesis vanishes for some 
O. As (3 > 0, this happens for some (real) positive v' , if and only if 



m) 



V COSfj) 



<0, 



(33) 



where <j> is the angle between r and v. If (|33[l holds, then the range of v' 
which a nonzero solution to (|32l) for A exists is 



for 



< v' < 



/?(fi) + 



v cost 



(34) 



(This is true for more-than-one dimensional space. If the space is one-dimensional, 
v' has only one acceptable value, as the parenthesis in 132|) has only two values 
at most one of them can be zero.) 

If (|33() holds, then there exists eigenvalues E for ft, with E > — v • v/4. 
Otherwise, all of the eigenvalues of ft are less than or equal to — v • v/4. The 
relaxation time of the system is 



1 



E n 



(35) 



where E max is the largest eigenvalue of ft. The largest value of E is either 
v ■ v/4, or the value obtained from i|29|) for the largest value of v' . So, 



4 
v-v ' 



-{min[2/3(n)+t> cos0]} 2 ' 



13(0) + ISpL 
13(0.) + ^f^ 



>0 
<0 



(36) 



In the first case, the system is in the fast dynamical phase, in which the relax- 
ation time does not depend on the boundary condition. In the second case, the 
system is in the slow dynamical phase, in which the relaxation time is larger 
and does depend on the boundary conditions. This is the dynamical phase 
transition. 



5 Concluding remarks 

It was seen that the so-called voting model defined on a one-dimensional lattice, 
has a natural analog on a multidimensional continuum. It was seen that there 
are two kinds of phase transition, a static one corresponding to a discontinuous 
change in the behavior of the stationary profile of the system, and a dynamical 
phase transition corresponding o the relaxation of the system toward its sta- 
tionary state. The static phase transition is controlled by the bulk reactions, 
while the dynamical phase transition is controlled by the bulk reactions and 
the boundary conditions both. This is analogous to what seen for the case of a 
one-dimensional lattice. 

There are, however, differences. In the multidimensional case, the static 
transition occurs when the direction of the drift velocity is changed. This can 
happen without being necessary that the drift velocity vanishes. In the one- 
dimensional case, however, the static phase transition occurs only when the 
drift velocity passes zero. The reason is that in one-dimension, the only way to 
change the direction of a vector smoothly, is that the vector vanishes at some 
point. 



The second difference concerns the dynamical phase transition, to be more 
precise, the largest eigenvalues of the operator h defined in l|25|) . In the one- 
dimensional case and in the slow phase, there is only one eigenvalue greater than 
the largest eigenvalue corresponding to the fast phase. In the multidimensional 
case, however, in the slow phase the spectrum of h contains a continuous region 
the lower bound of which is the largest eigenvalue of h in the fast phase. This 
means that in the one-dimensional case and in the slow phase, there is a largest 
relaxation time and a gap between this and the next largest relaxation time, 
while in the multidimensional case, there is no such gap. 
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